Impurity induced coherent current oscillations in one-dimensional conductors 
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We study theoretically the electronic transport through a single impurity in a repulsive Luttinger 
liquid (LL), and find that above a threshold voltage related to a strength of the impurity potential the 
DC current / is accompanied by coherent oscillations with frequency / = I/e. There is an analogy 
with Josephson junctions: the well-known regime of power-law I-V curves in the LL corresponds to 
damping of the Josephson current below the critical one, while the oscillatory regime in the LL can 
be compared with the Josephson oscillations above the critical current. 
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Basic electronic properties of three-dimensional (3D) 
solids are usually well described within Landau's Fermi- 
liquid picture where low-energy excitations are quasipar- 
ticles that in many respects behave like non-interacting 
electrons. This is not the case in ID systems where the 
usual Fermi liquid picture breaks down. In one dimension 
single-electron quasiparticles do not exist, the only low 
energy excitations are charge and spin collective modes. 
Such a state called the Luttinger liquid (LL) is an alter- 
native to Fermi liquid in ID (for a review see Ref. [l[). 
There are different realizations of the ID electronic sys- 
tems. The examples are semiconductor-based quantum 
wires in which dimensionality of the conduction elec- 
trons is reduced by dimensional quantization [21], metal- 
lic linear chains on Si surfaces Q, carbon nanotubes 
conducting polymers [H], and quantum Hall effect edge 
states ja]. There are also evidences that effects related 
to inter-electronic interaction can be taken into account 
in terms of LL in strongly anisotropic quasi-lD conduc- 
tors Q where the LL state can be stabilized by defects Q 
or by formation of the CDW gap induced by electron- 
phonon coupling [9<] . The transport properties of the LL 
are also very different from those of the Fermi liquid. 
In particular, isolated impurities form effectively large 
barriers in ID systems with repulsive inter-electronic in- 
teraction and strongly suppress the current which leads 
to a power-law dependence of conductivity 13, 11, 12] • 
This effect can be described in terms of tunnelling be- 
tween different minima of a periodic potential describing 
interaction of the electronic system with the impurity. 
The periodic potential is associated with Friedel oscilla- 
tions induced by impurity. The Hamiltonian describing 
interaction of the LL with the impurity in spinless LL is 
expressed in terms of the bosonic (plasmon) displacement 
field $(t,x) at the impurity position [lj 
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dxWi5(x) cos 2$>, 
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where Wj is related to 2fci?-component of the impurity 
potential. And the particle density operator reads [l[ 
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where the first term describes smooth variations of the 
particle density and the second one yields a rapidly oscil- 
lating part. Fluctuations of the field $ make the expec- 
tation value of the second term in the free LL equal to 
zero. However, fluctuations of $ are finite at the impurity 
which results in Friedel oscillations, i. e., in 2kp modula- 
tion of charge density decaying with distance as 
where K p is the LL parameter measuring the strength of 
the interaction: K p < 1 for repulsive and K p > 1 for at- 
tractive interaction, so the larger electronic repulsion the 
slower decay of the oscillations. As the current operator 
in the LL reads I — (e/ir)dt&, the current flow through 
the impurity implies an increase of $ with time, which 
corresponds to a shift of the Friedel oscillations. 

Power-law I-V curves are induced by tunneling be- 
tween minima of a washboard potential ([I]) slightly in- 
clined by an external bias, so that an increase of the 
phase by tt corresponds to a transition of one electron 
through the impurity [lj. This resembles Josephson junc- 
tions where similar fluctuations result in a finite voltage 
drop at a current below the Josephson critical current, 
I < I c (for a review see Ref. 13]). However, in supercon- 
ducting junctions, at I > I c the Josephson oscillations 
start, and this corresponds to an increase of the super- 
conducting phase difference with time in the washboard 
potential which is inclined to an amount exceeding the 
critical value. Below we show that a similar regime must 
occur in the LL with an impurity provided the applied 
voltage exceeds a threshold value corresponding to the 
slope at which the system can roll out from the min- 
imum of the washboard potential. Above the thresh- 
old the current is larger than the tunneling current in 
the sub-threshold regime of power-law I-V curves, and 
the current flow is accompanied by oscillations with the 
washboard frequency / = I/e. 

For brevity we consider the response of a spinless re- 
pulsive LL to an external DC voltage using Tomonaga- 
Luttinger (TL) model with short range interaction char- 
acterized by constant K p < 1. The short-range interac- 
tion describes gated quantum wires where the long-range 
part of the interaction is screened by 3D gate electrodes. 
At the end we will discuss essential modifications induced 
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by spin and the long-range Coulomb interaction. 

We start from the Hamiltonian that includes the stan- 
dard TL Hamiltonian, the impurity part ([T]), and the 
term with an external electric field, He = — J dx^E<t, 



H = dx 



1 



7T 2 K 2 



Hi + He 



(2) 



Commuting $ with the Hamiltonian we derive the equa- 
tion of motion for the Heisenberg operator 
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2W l sm2$ (t)6(x) -E , (3) 



where $o(t) = $(t,x = 0), and the operator in the left- 
hand side is the inverse free bosonic Green's function of 
density fluctuations. For the standard LL it reads 
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(v 2 c%-$)/(irv F ), 



where v = vpjK p is the velocity of plasmons. This op- 
erator does not contain damping. If one takes into ac- 
count coupling of electrons to a dissipative bosonic bath 
(to phonons or to density fluctuations in a metallic gate) 
then a finite damping appears, and the Fourier trans- 
formed retarded Green's function acquires the form [lit ] 
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As it was shown recently the damping suppresses soliton- 
like fluctuations and reduces power-law conductance to 
an exponential one [15| . 

Using Eq. §3§ with proper boundary conditions at the 
contacts one can express <fr(t,x) in terms of its value at 
the impurity site. We apply boundary conditions [16J 
derived for a wire adiabatically connected to ideal Fermi- 
liquid reservoirs at x = ±L/2 with voltage difference V. 
Then by means of the Fourier transformation we finally 
obtain the equation of motion for the operator &o(t) at 
the impurity position. 

d t $o(t) + - / dt 1 Z(t-t 1 )sm2$ (t 1 ) = -1,(5) 



Z(t) 



dio _ iut Wik p (l-ik p tan^) 



2tt 



\/uj 2 + iojv ~ j: 
qo = , K p = K p 



w + IV 



where / = e 2 EvF /irfiv is the time averaged macroscopic 
current in the channel. Eq. ([5]) is supplemented by rela- 
tion between time averaged values 



G(V-Vi), G- 1 = Go 1 +G-\ G„ 
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where Vi — 2Wi(sm2$o(t)) t is the DC component of 
the voltage drop at the impurity, Go = e 2 /h is the con- 
ductance quantum per spin orientation, and G v is the 



conductance related to damping in the clean part of the 
wire. 

As we want to concentrate on conduction through the 
impurity and not on the problem of contacts, we will 
consider the case of small enough damping, v <C ui, and 
a long channel, L ^S> l v = vjv. This allows to neglect 
reflections of current pulses generated by the impurity 
from the contacts and to substitute tan %p — > i. Then 
the kernel Z simplifies 



Z{t) = W i K p J ge" 
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which at v = gives Z(t) — * WiK p S(t). Note, however, 
that at any small but finite damping J Z(t)dt = 0, so the 
damping cannot be neglected at lj = 0. 

Remember that $o(0 m © is an operator, so it is 
not easy to solve this non-linear equation. So we extract 
the expectation value &o(t) — ($o), thus <lo = $o + 8$, 
where (<5$) = 0, and &o(t) satisfies the equation 

6 f°° 7T 

dt$o(t) + T / diiZ(t-ti)(cos25$(ti))sin2$o(*i) = -I- 



(8) 

To solve this equation one must calculate (cos2<5$) first. 
In calculation of this expectation value we will ignore 
the soliton-like fluctuations that are responsible for sub- 
threshold tunneling, but take into account Gaussian fluc- 
tuations that substantially reduce (cos2<5$). These fluc- 
tuations can be treated by means of the self-consistent 
harmonic approximation. 

But before treating the general case we discuss solu- 
tion in the simple limit of strong inter-electron interac- 
tion (K p — > 0) when Eq. ([5]) can be easily solved ana- 
lytically. In this limit fluctuations of the displacement 
field at the impurity are small, (cos 2S<t) — > 1, and <I>o(i) 
can be treated as c- number. When V < Vr = 2W% 
we find from Eqs. (J6j) and (jHJ) a stationary solution 
2<f>o = arcsin(VyVr) with zero current 1 = but with 
non-zero voltage drop over impurity. Note that we ob- 
tain that the current is zero in the sub-threshold region 
because we neglected solitonic fluctuations resulting in 
power-law I-V curves. 

At V > Vt the solution is oscillatory with fundamental 
frequency / = I/e 
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where In = GqVtK p . Eq. ([9]) determines the current at 
the impurity site, I(t,x = 0) = edt&o(t)/ir. Current at 
the clean part of the channel calculated from Eq. f3| is 
equal to I(t,x) = ed t &o(t — \ X \/ V )/ 7T at \x\ <C l v , and 
I(t, x) = I at large distances from the impurity. 
For the DC voltage drop at the impurity we obtain 



V t = V T 



(10) 
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Thus from ((6|), ([9]) and (fTO]) we see that the oscillatory 
regime starts at V > Vr- The current at the impurity 
consists of narrow pulses of height 2I at I -C Iq and 
transforms into the Ohmic current accompanied by har- 
monic oscillations of amplitude Iq at / 3> Iq. 

Now we consider the case of finite values of K p . Mean 
square fluctuations (<5<3E> 2 ) can be calculated from Keldysh 
Green's function D K = -i({6$(t),5$(t')} + ) at t = i! . 
This function can be expressed via retarded and ad- 
vanced Green's functions 

D R(A) = ±i6 [±( t _ t ')]([5$(t), 5$(t')]-) 

by relation D K (t,t') = D R {t,t')f(t') - f(t)D A (t,t'), 
where Fourier transform of / is related to the distribution 
function of bosonic excitations N(uj), f(ui) = 1 + 2N(lo). 
In the equilibrium state N(u>) is the Planck distribution. 
At low temperatures, smaller than all characteristic en- 
ergies of the system, one can neglect contribution of ther- 
mally excited excitations, then f(cu) — sign(w) and D K 
can be expressed via the retarded and advanced func- 
tions. Note that acting in this way we neglect the effect 
of non-equilibrium distribution of bosonic excitations. 

Now we derive equations of motion for the retarded 
and advanced Greens functions of fluctuations. This can 
be done in a standard way multiplying Eq. [3] by $(t', x') 
from the left and from the right, and combining them 
in order to obtain corresponding Green's function after 
averaging. Then using the Fourier transformation we ex- 
press Greens functions at the impurity site and get rid 
of the coordinate dependence of Green's function. After 
that by means of ([8]) we subtract expectation values and 
use the self-consistent harmonic approximation 

sin 25$ -► (cos25$)2<5$ = e" 2<54,2> 2(5$, 

and arrive, finally, to close equations of motion for D R ^ A ' . 
As we will need this equations at frequencies larger than 
the small damping constant we neglect v and employ 
Z(t) = W l K p 5{t -ti). Then the equation for D R reads 

d t D R (t,t>) + -W t K p C{t)D R {tX) = ^S(t,t>), 



we calculate the threshold voltage. In the sub-threshold 
region C does not depend on time. Then calculating 
integrals in (fTB")) and using definition of C (fTTj) we obtain 

(5$ 2 ) = ^£ In - (14) 

[ ' 2(1 - K p ) 2K p eWi cos 2$ ' 1 > 

where A ~ ppv is a large cut-off energy. Note that in 
accordance with our previous statement we found that 
fluctuations vanish at K p — ► 0. Substituting (fT4"|) in © 
we obtain equation for $o m the sub-threshold regime 



2Wi 



2K p eW l cos 2$ 
A 



Bin 2$o = V. (15) 



It has solution at V < Vt with 



V T = 2Wi 



' 2K 3 p /2 eW l \ 1 Kp 



V^K P - (16) 



C(t) = cos2$(ti)(cos2J$(ti)). 
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We see from (|16|) that the finite threshold voltage exists 
provided K p < 1, and this is in accordance with the 
condition that the impurity is a relevant perturbation in 
case of repulsive inter-electronic interaction. 

It is not simple to find (5<t 2 ) analytically in a general 
time-dependent case, but this can be done easily when 
the voltage slightly exceeds the threshold value and the 
DC part of the current is small, I <C Iq- Then $o in- 
creases with time in a step- like manner spending the most 
part of the oscillation period e/I near the value corre- 
sponding to maxima of C(t) and passes rapidly other 
values of C during short time interval ~ 1/ \JTTq that 
gives small contribution to the period. Therefore, in the 
most part of the period one can use (|14p with time depen- 
dent value <&o(t), and self-consistent solution of Eqs. (fT3"|) 
and ((SJ in this limit yields equations similar to ([9]) and 
(fPTJf with the value of the threshold voltage Vr (fl6|) . 

So far we considered the spinless LL. In the spinful 
LL impurities partly violate the spin-charge separation, 
and the impurity Hamiltonian contains the spin phase 
field [l[. This leads to modification of the results, in par- 
ticular, of the threshold voltage. For a spin-independent 
electronic repulsion Eq. (TT6|) is substituted by 



This equation and the similar equation for the advanced 
Green's function can be easily solved analytically. The 
solutions are 



itK, 



p0[±(t-t')}e T ^ mKp ^ C{tl)dtl . (12) 



D R(A) 



This gives us for (5$ 2 ) = \D K (t,t) 



{6 &) = f^ dh e -^K p ^ c{t-t 2)d t^ (13) 
2 Jo t\ 



Eq. (|T3|) must be solved self-consistently with {8J and 
© in order to find (<5$ 2 ) as function of cos2<I>o- First, 



1 + A' p 



2(1 -K p ). (17) 



Now we discuss the case of the long-range Coulomb 
repulsion. It can be taken into account in terms of mo- 
mentum dependence of K p (l7| 



K~ J = l+a*K {qd), 



4e 2 



where the McDonald function originates from the Fourier 
transformation of the Coulomb potential in a wire of di- 
ameter d, e is a background dielectric constant. In case 
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of the long-range interaction the impurity is a relevant 
perturbation at all values of the parameter a describing 
the strength of the repulsion. This is evident from the 
flow equation demonstrating that the impurity potential 
scales to infinity under renormalization 



dl 



1 + a 2 In • 
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The McDonald function is substituted here by its limiting 
logarithmic form. 

The momentum dependence of the interaction param- 
eter modifies equations and (jHJ) . With the logarithmic 
accuracy and for the case of negligible damping we find 
that Z(t) must be substituted for 



Z(t) = W t K p J g £ 



{uj + iu)(l + a 2 LnfH) 



(18) 



Though it is not simple to solve such equations analyti- 
cally we find that they do not result in important quali- 
tative difference from the case of the short-range interac- 
tion with constant K p . One of the most important dis- 
tinctions is that suppression of (008 2(5$) by fluctuations 
is smaller than in the TL model, and this results in the 
smaller effect of fluctuations on the threshold voltage. 
Calculating again (5$ 2 ) by means of Keldysh Green's 
function we find in the sub-threshold regime 



(5& 



du> 



2{e/h)WiC 



(19) 



Assuming a moderate strength of the Coulomb repulsion 
when a is of the order one, and taking into account that 
the argument of the logarithm in (TT5]) is large, we es- 
timate (<5<& 2 ) with the logarithmic accuracy. Then we 
obtain 



(6&) = -Jin 



edW.C 



Here we neglected the contribution from high-energy cut- 
off, assuming that In ^ ~ In is not a large value. 
This gives an estimate for the threshold field 



^ = 2^exp(-^ln — 



where, again, the argument of the square root in the ex- 
ponential function is given with the logarithmic accuracy. 

In conclusion, we found that above the threshold volt- 
age Vt the current through an impurity generates coher- 
ent oscillations with the fundamental frequency / = I/e. 
We considered the case of the DC applied voltage. If in 
addition there is also an AC component of the applied 



voltage with frequency fo then an analog of the Shapiro 
steps observed in Josephson junctions will appear on the 
I-V curves. In our problem these are the steps of a con- 
stant voltage, the fundamental step being located at the 
current value I = e/o- Characteristic frequencies of the 
oscillations induced by the DC voltage are determined 
by the strength of the impurity potential. In semicon- 
ducting quantum wires typical values of the impurity po- 
tential (say, from the shallow impurities) can be of the 
order of several meV, and depending on the strength of 
the electronic repulsion the frequency may fall into the 
gigahertz or terahertz frequency regions. Direct applica- 
tion of our results to real systems is limited by voltages 
smaller than distances to other electronic subbands. The 
results can be modified also due to different coupling of 
the ID system to 3D environment. 
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